Extraction of the Neutron Magnetic Form Factor from Quasi-elastic 3 He(e, e') at 

Q 2 = 0.1-0.6 (GeV/c) 2 



B. Anderson, 10 L. Auberbach, 18 T. Averett, 3 W. Bertozzi, 11 T. Black, n -0 J. Calarco, 22 L. Cardman^ 19 
G. D. Cates, 14 -0Z. W. Chai, 11 -© J. P. Chen, 19 Seonho Choi, 18 '© E. Chudakov, 19 S. Churchwell, 4 '© 
G. S. Corrado, 14 C. Crawford, 11 D. Dale, 21 A. Deur, 20 H P. Djawotho, 3 '! D. Dutta, n -0 J. M. Finn, 3 
H. Gao, n '0R. Gilman, 16 ' 19 A. V. Glamazdin, 9 C. Glashausser, 16 W. Glockle, 15 J. Golak, 8 J. Gomez, 19 
V. G. Gorbenko, 9 J.-O. Hansen, 19 '! F. W. Hersman 22 D. W. Higinbotham, 23 '! R. Holmes, 17 
C. R. Howell, 4 E. Hughes, 1 B. Humensky, 14 -0 S. Incerti, 18 'F|_ C. W. de Jager, 19 J. S. Jensen, 1 '! X. Jiang, 16 

C. E. Jones/'B M. Jones, 3 '! R. Kahl, 17 H. Kamada, 15 !^ A. Kievsky, 5 I. Kominis, 14 W. Korsch, 21 
K. Kramer, 3 -© G. Kumbartzki, 16 M. Kuss, 19 'B E. Lakuriqi, 18 '© M. Liang, 19 N. Liyanage, 19 '© J. LeRose, 19 
S. Malov, 16 D. J. Margaziotis, 2 J. W. Martin, 11 -© K. McCormick 13 -0 R. D. McKeown, 1 K. Mcllhany "'0 
Z.-E. Meziani, 18 R. Michaels, 19 G. W. Miller, 14 J. Mitchell, 19 '! S. Nanda, 19 E. Pace, 7 T. Pavlin, 1 '© 
G. G. Petratos, 10 R. I. Pomatsalyuk, 9 D. Pripstein, 1 '© D. Prout, 10 R. D. Ransome, 16 Y. Roblin, 20 '! M. Rvachev, 11 
A. Saha, 19 G. Salme, 6 M. Schnee, 18 J. Seely, 11 T. Shin, n >HK. Slifer, 18 -© P. A. Souder, 17 S. Strauch, 16 '© 
R. Suleiman, 10 '© M. Sutter, 11 '! B. Tipton, 11 Q L. Todor, 13 ^© M. Viviani, 5 B. Vlahovic, 12 ' 19 J. Watson, 10 
C. F. Williamson, 11 H. Witala, 8 B. Wojtsekhowski, 19 F. Xiong, n '0W. Xu, n >@ J. Yeh, 17 and P. Zolnierczuk 21 -© 

(Jefferson Lab E95-001 Collaboration) 

'California Institute of Technology, Pasadena, California 91125 
2 California State University, Los Angeles, Los Angeles, California 90032 
'College of William and Mary, Williamsburg, Virginia 23187 
* Duke University, Durham, North Carolina 27708 
5 Istituto Nazionale di Fisica Nucleare, Sezione di Pisa, 1-56100 Pisa, Italy 
6 Istituto Nazionale di Fisica Nucleare, Sezione di Roma, 1-00185 Roma, Italy 
' Universita di Roma "Tor Vergata" and Istituto Nazionale di Fisica Nucleare, Sezione di Roma II, 1-00133 Roma, Italy 
8 M. Smoluchowski Institute of Physics, Jagellonian University, PL-30059 Cracow, Poland 
9 Kharkov Institute of Physics and Technology, Kharkov 310108, Ukraine 

10 Kent State University, Kent, Ohio 44^4^ 
11 Massachusetts Institute of Technology, Cambridge, Massachusetts 02139 
ls North Carolina Central University, Durham, North Carolina 27707 
Old Dominion University, Norfolk, Virgina 23508 
1J> Princeton University, Princeton, New Jersey 08544 
15 Ruhr-Universitdt Bochum, D-44 780 Bochum, Germany 
16 Rutgers University, Piscataway, New Jersey 08855 
17 Syracuse University, Syracuse, New York 13244 
18 Temple University, Philadelphia, Pennsylvania 19122 
19 Thomas Jefferson National Accelerator Facility, Newport News, Virgina 23606 
20 Universite Blaise Pascal/IN2P3, F-63177 Aubiere, France 
21 University of Kentucky, Lexington, Kentucky 40506 
22 University of New Hampshire, Durham, New Hampshire 03824 
23 University of Virginia, Charlottesville, Virginia 22903 
(Dated: 20 November 2006) 

We have measured the transverse asymmetry At' in the quasi-elastic 3 He(e, e') process with high 
precision at Q 2 -values from 0.1 to 0.6 (GeV/c) 2 . The neutron magnetic form factor Gm was ex- 
tracted at Q 2 -values of 0.1 and 0.2 (GeV/c) 2 using a non-relativistic Faddeev calculation which 
includes both final-state interactions (FSI) and meson-exchange currents (MEC). Theoretical un- 
certainties due to the FSI and MEC effects were constrained with a precision measurement of the 
spin-dependent asymmetry in the threshold region of 3 He(e, e'). We also extracted the neutron 
magnetic form factor G M at Q 2 -values of 0.3 to 0.6 (GeV/c) 2 based on Plane Wave Impulse Ap- 
proximation calculations. 



PACS numbers: 13.40.Gp, 24.70.+S, 25.10.+S, 25.30.Fj 

I. INTRODUCTION 

The electromagnetic structure of the nucleon has long 
been a topic of fundamental interest in nuclear and parti- 
cle physics. First-order nucleon electromagnetic proper- 



ties are commonly parameterized in terms of elastic form 
factors 1]. At low values of four- momentum transfer 
squared, Q 2 , these functions have a simple interpreta- 
tion as the Fourier transforms of the nucleon charge and 
magnetization densities in the Breit frame. Their precise 
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experimental determination is important both for test- 
ing fundamental theories of hadron structure and for the 
analysis of other experiments in the field, such as parity 
violation measurements that are designed to probe 
the strangeness content of the nucleon. 

The proton form factors have been determined with 
good precision at low Q 2 using Rosenbluth separation of 
clastic electron-proton cross sections, and more recent! 



lently 
4,5]. 



at higher Q using a polarization transfer technique 
The neutron form factors are known less well because of 
the zero electric charge of the neutron, causing its electric 
form factor to be small, and experimental complications 
such as the lack of free neutron targets and difficulties 
associated with neutron detectors. 

Over the past two decades, with the advent of much im- 
proved experimental facilities, the precise measurement 
of both the neutron electric form factor, G^, and the 
magnetic form factor, G^ , has become a focus of activity. 
Until recently, most data on G^ had been deduced from 
elastic and quasi-elastic electron-deuteron scattering. In- 
clusive measurements of this type suffer from large the- 
oretical uncertainties due in part to the deuteron model 
employed and in part to corrections for final-state inter- 
actions (FSI) and meson-exchange currents (MEC). The 
sensitivity to nuclear structure is reduced by measuring 
the neutron in coincidence, 2 H(e, e'n) [6[, and, further, by 
taking the ratio of cross sections of 2 H(e, e'n) to 2 H(e, e'p) 
at quasi-elastic kinematics @, H, H, US, EH • Uncertainties 
of less than 2% in G^ have been achieved in the region 
Q 2 < 1 (GeV/c) 2 using the latter technique [IS [EH- De- 
spite this high precision, there is significant disagreement 



between the results of 
cent experiments 0, [lO 



I and those of the more re- 
llll ] of up to 10% in the absolute 



value of G^j. An explanation has been suggested in [12j, 
but the issue has remained contentious. 

To clarify the situation experimentally, additional data 
on Gjjf, preferably obtained using a complementary 
method, arc highly desirable. Inclusive quasi-elastic 
3 He(e, e0 scattering provides such an alternative ap- 
proach [l3j. In contrast to deuterium experiments, this 
technique employs a different target and relies on polar- 
ization degrees of freedom. It is thus subject to com- 
pletely different systematics. On the other hand, due 
to the more complex physics of the three-body system, 
the precise extraction of nucleon form factors from polar- 
ized 3 He measurements requires careful modeling of the 
nuclear structure and of the reaction mechanism. Re- 
cent advances in Faddeev calculations [3, EH 03 have 
brought theoretical uncertainties of 3 He models suffi- 
ciently under control to allow such studies in the non- 
relativistic kinematic regime. A precision comparable to 
that of the deuterium ratio experiments can be achieved 
using the polarized 3 He technique [l7l ]. 

The use of polarized 3 He targets was pioneered at MIT- 
Bates 0, HI M: H3 and Mainz [H. In [H, G n M was 
extracted for the first time from quasi-elastic inclusive 
scattering from polarized 3 He, although with a large sta- 
tistical uncertainty. 



In this paper, we report on the first precision measure- 
ment of the so-called transverse asymmetry At> , which is 
sensitive to G%j, in the inclusive reaction 3 He(e, e'). The 
results were obtained in Hall A at the Thomas Jefferson 
National Accelerator Facility (Jefferson Lab). Brief re- 
ports of these data have appeared previously [13, [H, [H| . 
This paper presents the data analysis and evaluation of 
model uncertainties in much more detail. In addition, the 
analysis has been slightly refined. The results presented 
here are final. 

The neutron magnetic form factor G^ was extracted 
12 = 0.1 to 0.6 (GeV/c) 2 in steps of 0.1 (GeV/c) 2 
|23j. In addition, high-precision asymmetry data in 
the 3 He breakup region were obtained at Q 2 -values of 
0.1 and 0.2 (GeV/c) 2 H3- The threshold data provide 
a stringent test of the above-mentioned Faddeev calcula- 
tions because they cover a kinematical region where the 
proper treatment of the reaction mechanism is particu- 
larly important. 

At the Q 2 = 0.1 and 0.2 (GeV/c) 2 kinematics, G n M 
was extracted using a state-of-the-art Faddeev calcula- 
tion [lq . At these low Q 2 , relativistic effects are small, 
and the non-relativistic Faddeev results have been shown 
to be in good agreement with a diverse set of few-bod 



)Ody 
24]. 



data, including our own He breakup threshold data 
On the other hand, the extraction of G^ f from our 3 He 
asymmetry data at higher values of Q 2 with the same pre- 
cision as that achieved at low Q 2 would require a more ad- 
vanced theory that includes both an accurate treatment 
of reaction mechanism effects (FSI and MEC) and proper 
relativistic corrections (and possibly other refinements, 
such as A-isobar excitations, presumed to be small at our 
kinematics). Unfortunately, such a comprehensive calcu- 
lation is not available at the present time, and efforts 
to extend the theory are only in the beginning stages. 
For example, full inclusion of FSI has been investigated 
for the two-body channel in [25J . The Hannover group 
has carried out a coupled-channel calculation of 3 He(e, e') 
that accounts for FSI and A-isobars [26j, unfortunately 
also with limited success at higher Q 2 . Nonetheless, we 
observe that the size of FSI and MEC corrections to in- 
clusive scattering data near the top of the quasi-elastic 
peak has been predicted to diminish sharply with increas- 
ing momentum transfer [53, [H, [H H3]. Hence, it ap- 
pears likely that the Plane Wave Impulse Approximation 
(PWIA), in which the knocked-out nucleon is described 
by a plane wave while the spectator pair is fully interact- 
ing, is reasonably accurate at the higher Q 2 -values of this 
experiment. A quantitative estimate of the Q 2 -behavior 
of deviations from the PWIA, in particular of the size 
of FSI corrections, could be obtained by performing a y- 
scaling analysis on the present 3 He asymmetry data |3l| . 
Such an analysis may be carried out in a future publica- 
tion. 

Taking the pragmatic point of view that the PWIA 
is currently the best available theory describing inclu- 
sive quasi-elastic scattering from polarized 3 He at Q 2 > 
0.3 (GeV/c) 2 , we have extracted G 1 ^ from our higher 
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FIG. 1: Spin-dependent inclusive electron scattering from a 
polarized target. The target spin angles, 6* and <j>* ar e defined 
with respect to the three-momentum transfer vector q. 



Q 2 data [23| using PWIA. While we do not attempt to 
go beyond the PWIA by computing corrections for the 
various effects omitted in this approximation, we provide 
estimates of the uncertainties of the results in consider- 
able detail. Despite the relatively large theoretical uncer- 
tainties in this approach, our results are in good agree- 
ment with the recent deuterium ratio measurements from 
Mainz [l(| [ll[ in the same <3 2 -region. 



II. THEORY 

A. Spin-dependent Inclusive Electron Scattering 

Figure [T] depicts inclusive scattering of longitudinally 
polarized electrons from a polarized nuclear target. The 
four-momentum of the electrons before and after the re- 
action is k = (E, k) and k' — (£?',k'), respectively. The 
four-momentum transfer to the target is q = k — k' — 
(u,q), with the usual definition Q 2 = — q 2 . 

The experiment measured the spin-dependent asym- 
metry A = (a + — er~)/(er + + <J~), where cr ± is the differ- 
ential cross section for quasi-elastic scattering of electrons 
with helicity h = ±1 from polarized 3 He. It can be ex- 
pressed in terms of nuclear response functions, R{Q 2 ,u>), 
and kinematic factors, v(Q 2 ,ui), as [32| 



A = 



cos9*v T 'Rt' +2sin#* cos ^vtl'Rtl 1 
vlRl + v T R T 



(1) 



where 0* and </>* are the polar and azimuthal angles of the 
target spin direction with respect to the three-momentum 
transfer vector, q, as shown in Figure Q] By choosing 9* 
— 0° or 6* = 90°, one can select the transverse asymme- 
try, At> , or the longitudinal-transverse asymmetry, Atl> ■ 

The nuclear response functions for inclusive quasi-el- 
astic scattering have been obtained through both PWIA 
and Faddeev calculations. These calculations will be dis- 
cussed briefly next. 



B. Plane Wave Impulse Approximation 

In the PWIA, it is assumed that a single nucleon within 
the target nucleus completely absorbs the momentum of 
the virtual photon and leaves the interaction region as a 
plane wave. The remaining two-nucleon subsystem still 
undergoes interaction. Exchange current effects are ig- 
nored. The target nucleus, in our case 3 He, however, 
is described by the solution of the Schrodinger equation 
with realistic nuclear forces. Relativistic effects are in- 
cluded by using relativistic energy conservation and a 
relativistic electron-nucleon cross section. 

The nuclear current tensor is calculated as the prod- 
uct of the nucleonic current tensor and the nuclear spec- 
tral function, which contains the nuclear structure in- 
formation (see for example [H, [H, HH, [36j]). The spin- 
independent part of the spectral function has the well- 
known interpretation as the probability of finding a nu- 
cleon of certain momentum and isospin in the target nu- 
cleus [l^l • The PWIA formalism available in the litera- 
ture is largely but not necessarily fully covariant. 

Expressions for the matrix elements of the nucleonic 
current tensor and the spin-de pen dent nuclear spectral 
function have been derived in [36| , The spectral func- 
tion can be computed numerically from the nuclear wave 
function, which in turn can be obtained from a model of 
the nucleon-nucleon (NN) potential. With the nucleonic 
current tensor and the nuclear spectral function at hand, 
expressions for the functions R(Q 2 ,ui), required in 
can be derived and evaluated numerically [36| , 

The PWIA results presented in the paper were cal- 
culated following [33| . The calculation was based on a 
3 He wave function derived from the Argonne AV18 NN 
potential [38| and used the Hohler nucleon form factor 
parameterization (39l |. 



C. Non-relativistic Faddeev Calculation 

In the Faddeev approach fl4|, the coordinate-space 
Schrodinger equation for three nucleons with two-nucleon 
interactions is decomposed into three separate equations 
[4q |. In momentum space, the three Faddeev equations 
can be written as three integral equations. The kernel 
in each equation involves only the interaction between 
one pair of the nucleons. Solutions are obtained numeri- 
cally. The Faddeev decomposition of the three-body (and 
four-body) problem has proven to be a very useful com- 
putational tool in studies of light nuclei. 

With regard to 3 He, the Faddeev formalism has been 
app lied to unpolarized pd and ppn electrodisintegration 
[l5l l4lj with full inclusion of all final-state rescatter- 
ing processes. This calculation was subsequently ex- 
tended to electrodisintegration of polarized 3 He [42[. A 
further extension was made by including proper treat- 
ment of meson-exch ang e currents [l6j according to the 
Riska prescription [421 , which relates NN forces and 
meson-exchange currents in a model-independent man- 
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FIG. 2: Meson-exch ane e current contributions included in the 
Faddeev calculation [TfJ . a) Couplings to a correlated nucleon 
pair; b) couplings to a tt or p in flight. 
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as a function of the neutron magnetic form factor, 
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(4) 



where \a\ ^> 1 and b > c at low Q 2 where the above 
assumptions hold. By comparing At> data with predic- 
tions for At' from a calculation, one can extract G^ f . 
The detailed procedure will be discussed in Section I VII 

For completeness we mention that, because \G P E \ ^> 
\G E \, the proton contribution to the transverse-longitu- 
dinal response, Rtl 1 , may be dominant despite the small 
effective proton polarization. Thus inclusive scattering 
from polarized 3 He is not a promising technique to mea- 
sure the neutron electric form factor, G E [2ll. |46|. 



III. EXPERIMENT 



ner through the continuity equation. In 16], only the 
dominant n- and p-like meson exchange terms shown in 
Figure were considered. The effect of A currents has 
also been studied and found to be small (see Section lVBj) . 

The derivation of the nuclear response functions in the 
Faddeev approach is described further in [13]. In this 
work, the resulting expressions were evaluated numeri- 
cally using the framework of [l6| for a large number of 
kinematical points corresponding to the acceptance re- 
gions covered by the experiment. The underlying 3 He 
wave function was obtained using the BonnB NN poten- 
tial Again, the Hohler parameterization [39j] was 
used to model the nucleon elastic form factors. The Fad- 
deev calculation does not include relativistic effects. 



D. Extraction of the Neutron Magnetic Form 
Factor 

Because the 3 He nuclear spin is carried mainly by the 
neutron, the spin-dependent response functions Rt> and 
Rtl' can be expected to contain a large if not dom- 
inant neutron contribution at quasi-elastic kinematics 
[l3| . Comparison of Equation (JTJ) with the corresponding 
expression for scattering from a free nucleon leads to the 
expectation (within PWIA) that 

R T , P n {G n M f + P P (G P M ) 2 (2) 
Rtu « P n G n M G n E + P p G p M G p E , (3) 

where P n and P p are the effective polarizations of the 
neutron and the protons, respectively, in 3 He. Because 
the proton spins largely cancel, we have \P p \ <C \P n \- 
Effecti ve p olarizations have been calculated e.g. in 
[HHE^. Since \G n M \ m \G P M \, the proton contribution 
to the transverse response Rt> is small, and hence Rt> is 
essentially proportional to (G^ f ) 2 . Based on these argu- 
ments, the asymmetry At' defined in (JT|) can be written 



A. Overview &z Kinematics 

The experiment, E95-001, was performed in Hall A at 
Jefferson Lab using a continuous-wave electron beam of 
15 /iA current and 70% longitudinal polarization, inci- 
dent on a high-pressure polarized 3 He gas target. The 
beam energies were 778 and 1727 MeV. 

Electrons scattered from the target were detected by 
two high-resolution spectrometers (HRS) positioned on 
the left and right-hand side of the beam line, respectively. 
Both spectrometers were configured for electron detec- 
tion and for independent operation (single-arm mode). 
The "electron spectrometer" on the left side of the beam 
performed the main physics measurement of inclusive 
3 He(e, e') scattering at six different quasi-elastic kine- 
matics. The second HRS, the "hadron spectrometer" to 
the right of the beam, detected 3 He(e, e) elastic scatter- 
ing and provided continuous high-precision monitoring of 
beam and target polarizations. The kinematic settings 
arc listed in Tabic fl] 



B. Polarized Electron Source & Beam Line 

The electron beam originated from a laser-driven 
"strained" GaAs source [43, [48j]. Polarized electrons 
were produced by illuminating a GaAs crystal in ultra- 
high vacuum with high-intensity circularly polarized laser 
light and removing electrons excited within the crystal 
by a strong external electric field. The polarization of 
the laser light was controlled electronically with the help 
of a Pockels cell. In this way, the electron beam helicity 
could be reversed rapidly (typically at 30 Hz), minimizing 
systematic errors in the measurement of spin-dependent 
asymmetries. To reduce systematic errors further, the 
overall sign of the beam helicity was reversed periodi- 
cally by inserting a half-wave plate into the injector laser 
light path. 
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The standard Hall A beam line instrumentation and 
beam raster [49| was employed. The beam energy was 
determined with an accuracy of better than 0.1% for all 
kinematics. 



C. Polarized 3 He Target 

The experiment employed an optically-pumped polar- 
ized 3 He gas target [49] of the spin-exchange type [50l |. 
The target cell of this system contained high-pressure 
(rs 10 atm) 3 He gas as well as admixtures of rubidium 
(to facilitate optical pumping) and nitrogen (to quench 
radiation trapping). While background from the rubid- 
ium was negligible, the nitrogen admixture contributed 
on the order of 10~ 2 to the total target number density, 
requiring a small dilution correction (see Section IIV CP . 

The target cell proper was a 40 cm long aluminum- 
silicate glass cell (p = 2.76 g/cm 3 ) with pa 1.2 mm thick 
walls and « 135 thick end windows. A second target 
cell, the so-called reference cell, was available for calibra- 
tion measurements. The reference cell had essentially the 
same dimensions as the target cell, except that it had no 
thin end windows but rather a uniform glass thickness 
throughout. Further details can be found in [5ll. I52T]. 

A typical 3 He nuclear polarization of 40% was a- 
chieved. The target spin direction was cither —62.5° ± 
0.5° or -243.6° ± 0.5° in the laboratory. (The differ- 
ence of the two angles was not exactly 180° because of 
a calibration inaccuracy.) The target spin was reversed 
regularly throughout the experiment to reduce system- 
atic errors from false asymmetries. 



Electron arm (quasi-elastic) 



Q' 2 


e e' e 


(GeV/c) 2 


(GeV) (GeV) (deg) 


0.1 


0.778 0.717 24.44 


0.193 


0.778 0.667 35.50 


0.3 


1.727 1.559 19.21 


0.4 


1.727 1.506 22.62 


0.5 


1.727 1.453 25.80 


0.6 


1.727 1.399 28.85 


Hadron arm (elastic) 




e e' e 


(GeV/c) 2 (GeV) (GeV) (deg) 


0.1 


0.778 0.760 23.73 


0.2 


1.727 1.691 15.04 



TABLE I: Kinematic settings for the quasi-elastic and elastic 
measurements. 



D. Spectrometers 

The two spectrometers were equipped with their stan- 
dard detector packages [49| consisting of a pair of Verti- 
cal Drift Chambers (VDCs) for tracking, two segmented 
scintillator planes to generate the trigger and provide 
time-of-flight information, and a CO2 gas Cherenkov de- 
tector for electron/pion separation. The HRSs had a us- 
able momentum acceptance of approximately 9%. For 
further pion rejection, a preshower and a total-absorption 
shower counter were employed in the electron-arm HRS, 
while the hadron-arm HRS was instrumented with two 
thin lead-glass shower counters. The geometric solid an- 
gle of each HRS was limited to 6.0 msr by a rectangular 
tungsten collimator. The central scattering angle was 
surveyed to better than 0.1 mrad. 

Trajectories of scattered particles were reconstructed 
using the VDC data and the standard optics model of 
the HRS [49J. The achieved momentum and scattering 
angle resolutions (a) were better than 0.05% and 2 mrad, 
respectively. The transverse (i.e. along the beam) posi- 
tion resolution at the target was approximately 2 mm. 

The pion rejection factor with the Cherenkov detec- 
tors alone was of order 100. Combining the Cherenkov 
and shower counters, a factor of over 1000 was achieved. 
Pion rejection was only a concern with the left-arm HRS, 
where pion production was not kinematically suppressed. 

IV. ANALYSIS 
A. Overview 

The experimental raw asymmetry was calculated as 



where N + and iV_ are the electron yields normalized by 
charge and electronic live time for positive and negative 
electron helicities, respectively. 

To extract the physics asymmetry, corrections had to 
be made for dilution, background, radiative effects, and 
bin centering. Sources of dilution were the finite beam 
and target polarizations, and scattering from the target 
walls and from the nitrogen gas in the target. Polarized 
background arose from the elastic radiative tail, which 
extended into the quasi-elastic region. Radiative correc- 
tions had to be applied to the raw quasi-elastic asymme- 
try. Bin centering corrections account for finite experi- 
mental acceptances. 

The normalized yields in (JSJ) can be written as 

N = N qe + N ert + N emp + N N * , (6) 

where N qe , N ert , N emp , and N N2 are the contributions 
of quasi-elastic scattering from 3 He (before radiative and 
bin centering corrections), the elastic radiative tail, tar- 
get wall ( "empty target" ) scattering, and scattering from 
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nitrogen in the target cell, respectively. Using ([6]), one 
can define dilution factors for each of the three back- 
ground contributions, 



R erap 
R N 2 
R ert 



jyemp 



J\fN 2 
jyert 



Nl e 



(7) 
(8) 
(9) 



and express the physics asymmetry as 



A phys = (I + Rert^ + R emp + R N 2 j 

-R erf A ert + AA qe + AA bm , 



j^exp 
PbPt 



(10) 



where PbPt is the product of beam and target polar- 
izations, A ert is the asymmetry of the clastic radiative 
tail, AA qe is the radiative correction to the quasi-elastic 
asymmetry, and AA bm , the bin centering correction. In 
Equation (TIT)]) , it is assumed that both the empty target 
and the N2 contributions have no asymmetry. During 
the analysis, the empty target and N2 false asymmetries 
were verified to be indeed consistent with zero. 

Among the various factors in (flT)]) , A exp , R emp , and 
R N2 could be determined directly from data, while R ert , 
A ert , AA qe and AA bm had to be determined from cal- 
culations or simulations. PbPt was monitored continu- 
ously during the experiment via elastic polarimetry and 
was determined as the ratio between the measured elas- 
tic asymmetry and the simulated elastic asymmetry, as 
described in Section [TV El 



B. Raw Asymmetries 

Raw asymmetries for both spectrometers were cal- 
culated according to Equation ([5]). The quasi-elastic 
data were analyzed in terms of electron energy loss, 
w = E — E', and grouped in bins of 10, 20, or 18.75 
MeV width, depending on Q 2 . The elastic data from the 
right-arm spectrometer were analyzed in terms of excita- 
tion energy, defined as 



E x = \ M 



1/T2 



2M{E-E') -4EE'sm 2 {0/2)-M, (11) 



where M is the mass of the 3 He nucleus and 9 the mea- 
sured electron scattering angle. The raw elastic asym- 
metry was obtained from the region —1 MeV < E x < 
+1 MeV. 

The angle between momentum transfer and target 
spin, 6* in Equation ([1]), varied between 0.2° and 10.0° 
depending on Q 2 . This resulted in an Rtl> contribution 
to the experimental asymmetry of less than 2%, as esti- 
mated by a PWIA calculation. The Rtl 1 contribution 
is included in the theoretical calculations that were used 
to extract G%[. Even though theoretical predictions of 



Rtl 1 are less accurate than those of Rt' (because of the 
uncertainty in G^), the uncertainty in our extracted 
due to Rtl 1 is negligible. 

Raw asymmetries obtained for the four different com- 
binations of of target spin orientation and overall beam 
helicity sign were compared to check for false asymme- 
tries. No statistically significant false signal was found. 
For the main physics analysis, data from the four po- 
larization configurations were combined to minimize the 
statistical uncertainty. 



C. Empty Target and Nitrogen Dilution Factors 

Because the target cell was sealed, background from 
the target cell wall could not be measured directly by 
emptying the target. In addition, the background rate 
from the nitrogen buffer gas in the target could not be 
easily calculated because the nitrogen partial pressure 
could only be determined approximately when the cell 
was filled. Therefore, it was necessary to determine both 
background yields in separate calibration runs with the 
reference cell. 

For each kinematics, quasi-elastic data were taken with 
the reference cell empty and filled with N2 at several pres- 
sure values. The reference cell nitrogen yield as a func- 
tion of nitrogen pressure was determined by subtracting 
the empty cell yield from the raw yields of the nitrogen 
runs. As the reference cell had physical dimensions very 
similar to those of the target cell, the reference cell ni- 
trogen spectra could be used as a direct measure of the 
target cell nitrogen yield, N N2 , provided that they were 
scaled to the nitrogen pressure inside the target cell. 

The nitrogen partial pressure in the 3 He target cell was 
determined as follows: As shown in Figure [31 the elastic 
nitrogen peak was clearly resolved in both the reference 
cell nitrogen spectrum (upper panel) and the spectrum 
from the 3 Hc target cell (lower panel), as measured with 
the right-arm spectrometer. As the nitrogen pressure cor- 
responding to the reference cell spectrum was known, the 
nitrogen pressure in the target cell could be determined 
by simple scaling. This procedure was only required for 
one kinematic setting since the nitrogen pressure was es- 
sentially constant throughout the experiment. The result 
waspAr 2 = 15.15±0.35 kPa. The variation of the N2 yield 
as a function of time was found to be within ±3%. We 
assigned an overall uncertainty of 5% to the measured 
nitrogen background yields. 

Obtaining the empty target yield (i.e. the yield due to 
scattering from the 3 He target cell walls) from the empty 
reference cell data was complicated by two factors: (1) 
the background yield from the cell walls was a function 
of beam position and the beam tune, and thus reference 
cell runs did not necessarily reflect the exact background 
conditions present during production data taking; and 
(2) the reference cell glass wall thickness and density were 
not equal to those of the target cell. 

Regarding (1), the variation of the empty target yields 
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obtained under nominally identical experimental condi- 
tions but at different points in time were compared and 
found to agree within ±15%. 

Regarding (2) we note (a) the target and reference cells 
were made of different types of glass, where the target 
cell glass density was about 9% larger than that of the 
reference cell; (b) the thickness of the reference cell side 
walls was found to be, on average, 2.5% thinner than that 
of the target cell, as determined by laser interferometry 
and (c) the target cell had very thin (135 /xm) end 
windows, while the corresponding reference cell end win- 
dows were about as thick (1.2 mm) as its side walls. End 
window contributions were minimized by using software 
cuts. The residual contribution of the thick end windows 
tends to compensate the effect of the thinner and less 
dense glass walls of the reference cell, although this is 
difficult to quantify. Hence, we assumed that the empty 
cell background yield of the 3 He target was identical to 
that of the empty reference cell (without the need for an 
explicit correction for the different cell properties) and 
assumed an overall systematic uncertainty of 25% in the 
empty yield, taking into consideration the statistical un- 
certainty, the time variation of the yield due to beam 
tune variations, and the differences in the cell properties. 

The empty target cell and the N 2 dilution factors 
(R emp and R N2 ) were determined by combining all empty 
target and nitrogen runs, respectively, at the same kine- 
matics. The nominator in Equations ([7]) and (JSJ) was cal- 
culated according to © as Nt e +N ert = N-N em P-N N2 . 
The time variation of the yields was included in the sys- 
tematic uncertainty of each contribution. The uncertain- 
ties are given in Table HTT1 

An ad-hoc upward correction of all the empty target di- 
lution factors by a factor of 2, which was used in our prior 
publications [It], [H, [57| , was dropped in this analysis as 
it had been motivated by an unphysical tail of appar- 
ently poorly reconstructed events seen in the right-arm 
spectrometer. Instead, a more conservative uncertainty 
was assigned to the empty target background subtraction 
at Q 2 = 0.1 and 0.2 (GeV/c) 2 , where the empty target 
background is largest. 



D. Monte Carlo Simulation 

A full Monte Carlo simulation was developed for this 
experiment [24j . which allowed averaging theoretical re- 
sults over the experimental acceptances and accounted 
for multiple scattering, ionization energy loss, external 
bremsstrahlung, and internal radiative corrections. 

To calculate the spin-dependent elastic and quasi- 
elastic radiative tails, internal radiation effects were mod- 
eled using the covariant formalism developed in [58[, 
generalized to the case of low-Q 2 quasi-elastic scatter- 
ing. This formalism accommodates polarization degrees 
of freedom. Standard, unpolarized radiative corrections 
[59[ | were applied to the elastic peak region. 



N 2 



800 




-4 -2 2 4 

H-Arm E x (MeV) 

3 He 



100 - 




-4 -2 2 4 



H-Arm E x (MeV) 



FIG. 3: Raw yields measured with the right-arm spectrometer 
in the region of the N2 elastic peak using the N2 reference cell 
(upper panel) and the 3 He target (lower panel) as a function 
of excitation energy, E x . The leftmost peak represents the 
N2 ground state, and the other peaks are related to excited 
states of N2. 
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E. Elastic Polarimetry 

The beam and target polarizations, P b and P tl were 
monitored continuously during the experiment using elas- 

I 



Here, the v% are kinematic factors, r = <3 2 /4M 3 2 Ho , and 
Ha = M 3 Ho(-^ 3 He/-^A r ) = —6.37. To allow direct com- 
parison with data, the Monte Carlo program described 
in Section [IV Dl was used to average Equation ([12]) over 
the experimental acceptance. We then obtained 

PbPt = Asira X fc/emp. (13) 
A el 

where A e ^ p and Af™ are the measured and simulated 
elastic asymmetry, respectively, and /n 2 and f em p are 
correction factors for the measured nitrogen and empty 
target cell dilution, respectively, for the elastic data sets. 

The data for A^f p are listed in Table HU Separate data 
are shown for each of the four possible spin and helicity 
configurations, which are largely consistent within their 
errors. For the evaluation of (|13p . the weighted average 
of the data for the four spin combinations was used. The 
dilution factors f^2 and f empty were obtained using the 
procedure described in Section HV CI 

No radiative corrections were applied to the elastic 
data since most radiative effects were included in the sim- 
ulation. Missing is the spin dependence of the Schwinger 
correction, which we deemed negligible. 

At the two beam energies, E = 0.778 and E = 1.727 
GeV, the overall relative systematic uncertainty in PbPt 
was 1.3% and 1.7%, respectively. In each case, the domi- 
nant contribution came from the uncertainty in the form 
factors F c and F m , followed by the contribution from the 
uncertainty in the target spin direction. 

The average P b P t so obtained was 0.208±0.001±0.004, 
where the errors are statistical and systematic, respec- 
tively. As a cross check, independent measurements of 
the polarizations were obtained using M0ller beam po- 
larimetry and NMR target polarimetry, yielding an over- 
all average value of P b P t = 0.215±0.013 [gj. The elastic 
polarimetry results were used for further analysis and av- 
eraged for each quasi-elastic kinematic setting separately 
(cf. Table [TTJ) to account for possible slow changes of 
the polarizations with time. The observed stability of 
the polarization data suggests that this procedure was 
adequate. 



tic polarimetry. As the 3 He elastic form factors, the 
charge form factor F c and the magnetic form factor F m , 
are known very well experimentally (r30l |. the 3 He elastic 
asymmetry can be calculated as [321 ] 



(12) 

I 

V. ASYMMETRY RESULTS 
A. Quasi-elastic Transverse Asymmetry A T / 

Results for the quasi-elastic transverse asymmetry At> 
at the six measured Q 2 -points are shown in Figure[4j Nu- 
merical values can be found in 57]. The errors on the 
data are statistical only, while the systematic uncertainty 
is shown as an error band at the bottom of each panel. 
A detailed breakdown of the systematic uncertainties is 
presented in Table IIIII The experimental data were cor- 
rected for radiative effects, background, and dilution, as 
described in detail in the previous section. 

Also shown in Figure Q] are the results of several cal- 
culations. Dashed lines represent the PWIA calculation 
[33| . The dash-dotted and solid curves at the two kine- 
matics with lowest Q 2 represent, respectively, Faddeev 
results with inclusion FSI only [15j and with inclusion 
of both FSI and MEC corrections[l6] • Calculation [l]| 
will be referred to as the "full Faddeev calculation" in 
the following. All theory results were averaged over the 
spectrometer acceptances using the Monte Carlo simu- 
lation described in Section HVDl Further details on the 
calculations are given in Section HT1 

One observes excellent agreement of the data with the 
full Faddeev calculation over the entire w-range at Q 2 = 
0.1 and 0.2 (GeV/c) 2 , while PWIA describes the data 
well at the higher Q 2 , in particular in the region around 
the quasi-elastic peak (near the center of the w-range in 
each panel). 

B. Asymmetry in the Threshold Region 

The asymmetries measured in the region around the 
two- and three-body breakup thresholds (5.5 and 7.7 
MeV, respectively) are shown in Figure [5] These results 
provide a sensitive test of the quality of the Faddeev cal- 
culations. 

The threshold asymmetry data were taken with the 
hadron-arm spectrometer as a by-product of the elastic 
polarimetry and were analyzed in the same manner as 
the quasi-elastic asymmetries. The kinematics are given 
in the lower panel of Table HI The Q 2 -values of 0.1 and 
0.2 (GeV/c) 2 in Figure [5] correspond to the momentum 



-2tv T ' cos6*h 2 a F^ + 2 v /2t(1 + t)v tl > sinfl* cos </)*n A ZF m F c 
(1 + t)v l Z^F^ + 2tv t ^ a F 7 1 
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ri 2 






\K7(%)\ 




(GeV/c) 


(GeV/c) 


—62.5° / in 


— 62.0 / out 


— 24o.6 / m 


— 243.6 / out 


0.1 


0.1 


1.333 ± 0.027 


1.043 ± 0.027 


1.067 ±0.02 


1.208 ± 0.030 


0.193 


0.1 


1.078 ± 0.037 


1.177 ± 0.027 


1.190 ±0.021 


1.102 ±0.023 


0.3 


0.2 


1.251 ±0.096 


1.222 ± 0.048 


1.107 ±0.067 


1.206 ± 0.075 


0.4 


0.2 


1.181 ±0.055 


1.314 ±0.061 


1.168 ±0.06 


1.258 ±0.057 


0.5 


0.2 


1.265 ± 0.042 


1.307 ±0.039 


1.200 ± 0.045 


1.184 ±0.041 


0.6 


0.2 


1.258 ±0.049 


1.301 ± 0.047 


1.110 ±0.05 


1.096 ±0.05 



TABLE II: The measured elastic asymmetries |^f p | for the six quasi-elastic kinematic settings. Q^ e and Q\i are the momentum 
transfers of the quasi-elastic and elastic measurements, respectively. The four columns of results correspond to the four 
combinations of the signs of the target spin and beam helicity. The column headings indicate the laboratory target spin angle 
and the position of the accelerator injector half-wave plate. The uncertainties are statistical. 
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FIG. 4: (Color online.) Quasi-elastic At 1 asymmetry results vs. the energy transfer ui. Errors on the data points are statistical. 
The systematic uncertainty is shown as an error band at the bottom of each panel. 



10 




PWIA (a) 
FSI only 

FSI+MEC(no A current) 
FSI+MEC(with A current) 
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PWIA (b) 
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FSI+MEC(no A current) 
FSI+MEC(with A current) 
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FIG. 5: (Color online.) The experimental asymmetry in the region of the 3 He breakup threshold together with theoretical 
calculations for (a) Q 2 = 0.1 (GeV/c) 2 and (b) Q 2 — 0.2 (GeV/c) 2 . The calculations differ only in the description of the 
reaction mechanism. 
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transfer at the elastic peak. The data are plotted as a 
function of the excitation energy, E x , defined in Equa- 
tion (fTT|) . Horizontal errors represent the uncertainty in 
determining E x , which was dominated by the uncertainty 
in the beam energy. The vertical errors arc the statistical 
and systematic errors added in quadrature. Tables of the 
data and uncertainties can be found in [62j. 

Figure [5] also shows various theoretical results. Dot- 
dashed lines depict those of the PWIA calculations 1331] . 
while those of Faddeev calculations with FSI only [15| 
appear as dashed lines. The "full Faddeev calculation", 
which includes both FSI and MEC, but not the A- isobar 
current, [HI is represented by dotted lines. The solid 
lines were obtained with the full calculation after includ- 
ing A-isobar currents. These calculations employed the 
AV18 NN interaction potential. Results obtained with 
the BonnB potential were found to be only slightly differ- 
ent from the AV18 results and in even better agreement 
with the data [24[ • 

As can be seen, the agreement between PWIA calcula- 
tions and the data is poor at both kinematics, which con- 
firms the expectation that FSI and MEC corrections are 
essential in this region. Indeed, the inclusion of FSI im- 
proves the agreement significantly, and good agreement 
is achieved once MEC are added. It has been shown that 
substantial MEC are needed to describe the measured 
clastic electromagnetic form factors of three-nucleon sys- 
tems [gU . The corresponding physics should extend into 
the low-w region of inelastic scattering as well. 

The good agreement between the full calculation and 
the data at Q 2 = 0.1 (GeV/c) 2 suggests that FSI and 
MEC are properly treated in the full calculation. The 
insensitivity of the results to the addition of A-isobar 
currents implies a weak model dependence of the MEC 
corrections. The small, systematic discrepancy at Q 2 = 
0.2 (GeV/c) 2 may indicate that some Q 2 -dependent ef- 
fects, such as relativistic and three-nucleon force effects, 
become important already at this momentum transfer. 



VI. EXTRACTION OF THE NEUTRON 
MAGNETIC FORM FACTOR 

The neutron magnetic form factor, G 7 },, can be ex- 
tracted from the measured 3 He quasi-elastic transverse 
asymmetry At* if a calculation is available that predicts 
At i as a function of G^ f . If we assume, following Equa- 
tion that the asymmetry is a function of (G^) 2 , we 
can expand At> around a reference value, Go, 



solved for Gl,, , assuming the second-order term is small: 



A T >(G n M 2 ) = A TI {Gl) + 



dAr, 



d(G 



+ 0((G^ 2 -G 2 ) 2 ). 



~(Go) x (G^,/ 2 - Gq) 
(14) 



For ease of notation, we normalize all G%j values to a 
convenient reference scale (the Holder parameterization 
[39l | in this case) so that Go = 1. Equation IT4]) can be 



Gn 
M — 



A T ,{G1 2 )~A T ,{1) 
dA T ./d{Gl 2 m ' 



(15) 



Here, At'{G%j ) is the measured asymmetry. The pre- 
dicted asymmetry, At>(1), and the sensitivity factor, 
dA T > /d(G 7 lj 2 )(l), are the output of the model calcula- 
tion. The latter two parameters were determined using 
the full Faddeev calculation [l6[ for the lowest two Q 2 
points of this experiment, and the PWIA calculation [33| 
for the remaining four Q 2 . Results were averaged over the 
experimental acceptance using the Monte Carlo program 
described in Section IIVDI At each kinematical point, 
asymmetries were generated for several w-bins around 
the quasi-elastic peak. Within each bin, G^ was varied 
around the reference value Go by adding a constant to 
the functional form G r l I (Q 2 ) given by the Hohler model. 

G%r (Q 2 ) was extracted for each w-bin via (fT5|) . A dif- 
ferent functional form, a general second-order expansion 
of At , (G'1 [ ) 1 was also tried. The differences between the 
form factors extracted via these two methods was found 
to be negligible (< 0.1%) for all kinematics [13] 

The final G^ results were obtained by taking the 
weighted average of the G^ values from the w-bins clos- 
est to the quasi-elastic peak. The w-region used for the 
extraction of G 7 ^ covered a width of 30 MeV at Q 2 = 0.1 



and 0.2, 60 MeV at Q 2 = 0.3, 40 MeV at Q 2 
0.5, and 56.25 MeV at Q 2 = 0.6 (GeV/c) 2 . 



0.4 and 



The extraction procedure gives rise to a systematic er- 
ror due to the uncertainty in the experimental determi- 
nation of the energy transfer uj (±3 MeV). The uncer- 
tainty in lu results in an uncertainty as to the w-region 
over which to integrate the theoretical calculation used 
for the extraction of GJf . A shift of bin boundaries gen- 
erally translates into a different average value of At> for 
the bin and hence a different extracted G^ value. 

Furthermore, as can be seen in Figure UJ the theo- 
retical calculations, especially PWIA, match the data 
best in the immediate vicinity of the quasi-elastic peak 
where corrections to the plane-wave picture are small- 
est, whereas deviations may occur off the peak. This can 
introduce an artificial w-dependence into the extracted 
G 7 Ij which goes beyond the effect of the kinematical vari- 
ation of Q 2 with u>. For this effect to be minimized, 
the bins used for the G^ extraction should be centered 
around the quasi-elastic peak, assuming that deviations 
are distributed roughly symmetrically. The experimental 
uncertainty in ui may cause improper centering, resulting 
in a bias in extracting G^. The calculated uncertainties 
in G^j resulting from the uncertainty in oj can be found 
in Table M 
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Source 


SA T > 


AM 


%) 


( 


2 2 < 0.2 


Q 


2 > 0.3 


PtP b 


1.3 




1.7 


Empty target subtraction 


1.0 




0.25 


N2 background subtraction 




0.3 




QE radiative correction 




0.3 




Elastic radiative tail 




0.3 




Spectrometer acceptance 




0.5 




HC scintillator efficiency 




0.1 




HC wire chamber efficiency 




0.1 




HC computer deadtime 




n 1 
U. 1 




HC beam current shift 




0.1 




HC beam motion 




0.1 




Pion contamination 




0.1 




Total 


1.8 




1.9 



TABLE III: Estimated systematic uncertainties of the quasi- 
elastic A T i asymmetry measurements. "HC" denotes "helici- 
ty-correlated" . The two columns of uncertainties correspond 
to the quasi-elastic measurements at lower Q 2 (0.1 and 0.2) 
and higher Q 2 (0.3 — 0.6), respectively. Values in the center 
of both columns are common to all kinematics. 
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The effect of different NN potential models on the pre- 
dicted asymmetry At< was studied by carrying out the 
full Faddeev calculation with the Argonne AV18 and the 
Bonn B NN potentials at several representative kinemat- 
ics. In a similar manner, to estimate the uncertainty due 
to the elastic nucleon form factors other than G^, Fad- 
deev calculations were performed in which these quan- 
tities were varied individually by their published exper- 
imental uncertainties. The resulting uncertainty in G 1 ^ 
from these sources, when combined in quadrature, is less 
than 1% for all kinematics (cf. Table fV*|) . 

B. Relativistic Effects 

Since the full Faddeev calculation is non-relativistic, 
it was particularly important to estimate quantitatively 
the size of relativistic corrections. An approximate esti- 
mate can be obtained within PWIA, which is theoreti- 
cally well understood. Standard PWIA calculations take 
most relativistic effects into account (cf. Section Ul Bp . It 
is straightforward to modify the relativistic parts of the 
PWIA formalism to reflect the non-relativistic approxi- 
mations made in the Faddeev formalism. The differences 
between the results of such a modified, non-relativistic 
PWIA calculation and the standard relativistic PWIA 
results provide an estimate of the error in the Faddeev 
results due to relativistic effects. 

To this end, we modified three parts of the standard 



FIG. 6: (Color online.) Relativistic effects in A T i . The 
solid line is the standard, relativistic PWIA calculation [33l ]. 
The dot-dashed curve is the non-relativistic PWIA calculation 
that we developed, and the dashed curve is the non-relativistic 
PWIA calculation with heuristic relativistic corrections ap- 
plied (see text). 



PWIA formalism: approximations were made to the rel- 
ativistic kinematics, the phase space and the integral 
ranges of the Fermi momentum and the missing mass of 
the many-fold integration of the 3 He(e, e') cross-section 
were changed according to the non-relativistic kinemat- 
ics, and the relativistic hadronic current was trans- 
lated into an approximate, non-relativistic form (r34l. [6o|. 
Among the three modifications, the change of the kine- 
matics was found to dominate [66j |. 

With the PWIA results at hand, we developed a heuris- 
tic "recipe" [66| to allow an approximate correction of the 
Faddeev results for relativistic effects. The "recipe" could 
be readily applied to existing Faddeev results without the 
need for recomputation. 

Results of these studies are shown in Figure O The 
three curves represent the original relativistic PWIA 
results (solid line), non-relativistic PWIA results ob- 
tained using the modifications described above (dot- 
dashed line) , and non-relativistic PWIA results corrected 
for relativistic effects though the "recipe" (dashed line). 
As can be seen, the heuristic correction works well up to 
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FIG. 7: (Color online.) FSI effect study. The dashed curve 
represents the standard (relativistic) PWIA calculation, the 
solid curve is the (non-relativistic) Faddeev calculation with 
FSI effects only, and the dot-dashed curve depicts the same 
Faddeev calculation but with relativistic corrections applied. 
Comparing the dashed and dot-dashed curves, one can esti- 
mate the effect of FSI in A T i . 

about Q 2 = 0.4 (GeV/c) 2 . 

The acceptance-averaged difference between the rela- 
tivistic and non-relativistic PWIA results at Q 2 = 0.1 
and 0.2 (GeV/c) 2 was taken as the model uncertainty of 
the Faddeev results due to relativity. 

C. FSI & MEC 

To estimate FSI contributions to At>, we carried out 
the Faddeev calculation up to Q 2 = 0.4 (GeV/c) 2 with 
the inclusion of FSI effects only. (Already at Q 2 = 
0.3 (GeV/c) 2 , the 3N center-of-mass energy is above 
the pion production threshold, and therefore the non- 
relativistic framework is no longer valid.) Next, we ap- 
plied relativistic corrections to the Faddeev results using 
the ad- hoc prescription developed in Section IVlI Bl The 
"relativistic FSI" results so obtained were compared to 
the results of the standard (relativistic) PWIA calcula- 
tion, as illustrated in Figure [7] The difference between 
the two calculations in the region around the quasi-elastic 
peak is a measure of the FSI effects at each Q 2 point. For 



FIG. 8: (Color online.) MEC effect study. Comparing the full 
calculation (solid curve) with the calculation with FSI effects 
only (dashed curve), one can estimate contributions to A T / 
from MEC effects. 



the two highest Q 2 values, we extrapolated the FSI data 
using a purely empirical fit to the lower Q 2 values, as 
shown in Figure ^a). As expected [H, S<| , FSI ef- 
fects decrease significantly as Q 2 increases. 

In a similar manner, we can estimate the size of MEC 
effects by comparing the Faddeev results with inclusion of 
FSI only, obtained in the FSI study above, to those of the 
full Faddeev calculation. Results are shown in Figure 
and differences between the two calculations are plotted 
as solid triangles in Figure[5Jb). As with FSI, we observe 
a sharp decrease of MEC corrections with increasing Q 2 . 

It is interesting to compare our results for the size of 
MEC corrections with those obtained from theoretical 
studies of quasi-elastic inclusive scattering from polar- 
ized deuterium, d(e,e') [30| . The deuterium results are 
shown in Figure^b) as solid squares. As can be seen, the 
data are similar to those calculated for the correspond- 
ing 3 He reaction. Assuming a similar underlying physical 
mechanism, we use the MEC data from deuterium to es- 
timate the size of MEC corrections to the 3 He(e, e') data 
at the highest two Q 2 values of our data set. 

These studies provide rough estimates of the expected 
magnitudes of the respective effects. They are not reli- 
able enough to be used to correct the PWIA results for 
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0.6 

Q 2 (GeV/c) 2 



FIG. 9: Estimated magnitude of (a) FSI effects and (b) MEC 
effects in A T , as a function of Q 2 . In the lower panel (b), 
the solid triangles represent results obtained in our study of 
3 He(e, e'), while the solid squares depict predictions from a 
calculation of d(e, e') obtained in [3p| ]. The curves are empir- 
ical fits to the data. 



FSI and MEC contributions. Consequently, we use the 
numbers obtained above as estimates of the model un- 
certainties inherent in the PWIA. We take the numbers 
as the la values of the uncertainties, which we assume to 
be symmetric. The resulting model uncertainties in 
are detailed in Table IVl and are propagated into the final 



G^f errors given in Table [TV 
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1.26 
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2.8 
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0.9577 


1.35 
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0.9689 


1.35 


1.25 
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2.8 


0.6 


0.9939 


1.55 


1.38 
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D. Off-Shell Effects 

Off-shell corrections to the single nucleon current, in- 
cluding the part of the current that describes polarization 
degrees of freedom [63, are purely relativistic in nature. 
While the PWIA calculation used here includes off-shell 
effects, they are ignored in the Faddeev calculation. 

We estimated the magnitude of required off-shell cor- 
rections to the Faddeev results by comparing results of 
a modified version of the PWIA formalism that treats 
nucleons as on-shell 66] to those of the standard PWIA. 

In addition, theoretical uncertainties due to different 
possible off-shell prescriptions were estimated using the 
difference of PWIA results obtained with the deForest 
CC1 and CC2 forms (The standard PWIA calcula- 
tion employs CC1.) While this number represents a min- 
imum uncertainty, as various other off-shell prescriptions 
are equally permissible (6?| , PWIA calculations using the 
CC1 form have been found to agree better with exper- 
imental data of unpolarized 3 He(e, e') scattering than 
those using other prescriptions 33]. This suggests the 
use of the CC1 prescription as a reference in the polar- 
ized case as well. 

Results are given in Table [V] Interestingly, off-shell ef- 
fects dominate the model uncertainty in G^ at the lowest 
two Q 2 -values. 



VIII. FORM FACTOR RESULTS AND 
DISCUSSION 



Numerical values for G 1 ^ extracted in this work are 



TABLE IV: Results for G%[ as a ratio to the dipole form 
factor, Gd, and uncertainties obtained in the present exper- 
iment. The data have changed slightly from our previously 
published numbers fvH . [23j] due to differences in the analysis. 



given in Table IIVI (in units of the empirical dipole pa- 
rameterization, Gd = (1 + Q 2 /0.71) -2 ) and shown in 
Figure [10] along with the existing world data set pub- 
lished since 1990 [E 0, 1, ©, M, El, HI . The error bars 
represent the total uncertainties reported by the respec- 
tive experiments, including model uncertainties. 

The results appear to be largely consistent, with the 
exception of the early 2 H(e, e'n) data from Bates @ and 
the first 2 H(e, e'n)/ 2 H(e, e'p) ratio measurement from 
Bonn 0, Q. The discrepancy of the data of these two 
experiments with the rest of the world data has been at- 
tributed to incomplete corrections for neutrons that miss 
the neutron detector The data of the Bonn exper- 
iment [3] were re-analyzed subsequently Q , resulting in 
a downward correction of the G 7 ^ data. Figure [TOl shows 
the re-analyzed data. We note the satisfactory agreement 
between the more recent, high-precision deuterium ratio 
measurements 0, OiJ [HI and the data from this work. 
The agreement is well within the total uncertainties of 
the experiments, except at Q 2 = 0.5 and 0.6 (GeV/c) 2 , 
where the 3 He results are low by about 5%. 

Also shown in Figure [TU] are several theoretical re- 
sults: a recent dispersion-theoretical fit by Hammer and 
Meifiner [691 (solid curve), a chiral soliton model by 
Holzwarth [70] (dotted curve), a relativistic baryon chiral 
perturbation theory calculation by Kubis and Meifiner 
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Q 2 (GeV/c) 2 



FIG. 10: The world's G n M data since 1990. Data points represent the results of the Bonn (□), MIT-Bates @,|2(| (0,O), 

NIKHEF/PSI || ("&) and the Mainz/PSI [ToL [Tl| (A,x) experiments as well as those of the present measurement (•), where 
the error bars are the total uncertainties reported. Also shown are the results of various model calculations: Hammer and 
Meifiner [69l| (s olid curve), Holzwarth [7(J (dotted), Kubis and Meifiner p7l| ] (long-dashed), Lomon [7^1 (dashed-dotted) , and de 
Melo et al. [73] (short-dashed). 
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0.95 0.15 0.38 0.60 0.38 


1.25 


0.40 


0.5 
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0.5 
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1.0 


0.5 2.1 


0.6 


0.95 0.10 0.32 0.64 0.69 


1.38 


0.40 


0.5 
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TABLE V: Estimated uncertainties of the extracted form factor G 7 ^. Systematic uncertainties include contributions from the 
asymmetry measurement (A T r, see table HIT)) . the energy transfer determination (u), and the other nucleon form factors (G E , 
G P M , and G E ). Theoretical (model) uncertainties include contributions from the NN potential model, off-shell effects, final-state 
interactions (FSI), meson-exchange currents (MEC), three-body forces (3BF), Coulomb corrections, and relativistic effects. In 
the totals, the uncertainties have been added in quadrature, ignoring any possible correlations between the contributions, which 
may very well exist, especially for the model uncertainties. Thus, the numbers should be taken with appropriate caution. 



[7l| (long-dashed curve), a vector meson dominance 
(VMD) fit by Lomon [72j (dashed-dotted curve), and 
a recent light-front quark model by de Melo et al. (73j 
(short-dashed curve). It should be noted that all of these 
models contain one or more free parameters that have 
been fitted to existing data. 

As can be seen, the dispersion-theoretical fit [69| and 
the VMD fit [zi| agree best with the data at Q 2 > 0.3 
(GeV/c) 2 , while the chiral perturbation theory (ChPT) 
results 7l| and the chiral soliton model [7C} match the 
data better at lower Q 2 . The ChPT model is expected to 
be good only up to Q 2 « 0.3 (GeV/c) 2 , but clearly works 
very well in its reg ion of validity. The light-front model of 
de Melo et al. [73| arguably shows the best overall agree- 
ment. The models by Holzwarth [7(|, Lomon [72|], and 



de Melo et al. [73| also describe the proton form factor 
ratio G P E /G P M and other elastic nucleon form factors well 
in this Q 2 -region. 



IX. CONCLUSIONS 

In conclusion, we have determined the neutron mag- 
netic form factor G^ from quasi-elastic 3 He(e, e') data. 
At Q 2 of 0.1 and 0.2 (GeV/c) 2 , we used a state-of-the- 
art Faddeev calculation that includes FSI and MEC, and 
PWIA at four additional points between Q 2 = 0.3 and 
0.6 (GeV/c) 2 . The results agree within the total un- 
certainties with those obtained by several recent mea- 
surements on deuterium, except at Q 2 = 0.5 and 0.6 
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(GeV/c) 2 where the 3 He results are slightly low. A con- 
sistent picture of the behavior of in this <3 2 -region 
is beginning to emerge, although further precision mea- 
surements as well as improved model calculations, such 
as the extension of the Faddeev formalism to higher Q 2 
in the case of polarized 3 He, remain highly desirable. 

In addition, we have measured At> in the two- and 
three-body breakup threshold region at Q 2 of 0.1 and 
0.2 (GeV/c) 2 where the sensitivity to FSI and MEC ef- 
fects is particularly high. The results agree well with 
the predictions of the Faddeev model, especially at Q 2 = 
0.1 (GeV/c) 2 , confirming the validity of the treatment of 
FSI and MEC effects in this formalism. 
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